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Abstract 

The ground-state properties of quartet quantum-dot arrays are studied using the numerical 
renormahzation group (NRG) method with a four-site Hubbard model connected to two non- 
interacting leads. Specifically, we calculate the conductance and local charge in the dots from the 
many-body phase shifts, which can be deduced from the fixed-point eigenvalues of NRG. As a 
function of the on-site energy which corresponds to the gate voltage, the conductance shows 
alternatively wide peak and valley. Simultaneously, the total number of electrons N^i in the four 
dots shows a quantized stair case behavior due to a large Coulomb interaction U . The conductance 
plateaus of the Unitary limit emerging for odd N^i are caused by the Kondo effect. The valleys of 
the conductance emerge for even Nf,\, and their width becomes substantially large at half- filling. It 
can be regarded as a kind of the Mott-Hubbard insulating behavior manifesting in a small system. 
These structures of the plateaus and valleys become weak for large values of the hybridization 
strength F between the chain and leads. We also discuss the parallel conductance for the array 
connected to four leads. 

PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 
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I. INTRODUCTION 



The Kondo effect in quantum dots is a subject of current interest, and early theoreti- 
cal predictions have already been confirmed by the experiments. j3|, y| Recently, the 
interplay of various effects such as Aharanov-Bohm, Fano, Josephson, and Kondo effects 
have been studied extensively, and interesting phenomena caused by 
inter-electron interactions have been expected to be seen at low temperatures. 

For studying the low-temperature transport of interacting electron systems, careful cal- 
culations are required. So far, some numerical approaches to conductance through small in- 
teracting systems have been examined by several groups. For instance, a Fermi-liquid based 
method, 3| density matrix renormalization method, 14 1 and functional renormaliza- 
tion method, 3 Q| seem to be complement each other to give reliable information valid 
at low temperatures. The Wilson numerical renormalization group (NRG) method, which 



successi vely eliminates higher-energy states to obtain accurately the low-lying energy states, 
[l^! 12, 3 has also been applied successfully to single and double quantum dots. P, 13 
Recently, we have provided an explicit prescription to deduce the many-body phase shifts 
from the fixed-point energy levels of NRG. It is applicable to a wide class of the quantum-dot 
systems connected to noninteracting leads. With this method, we examined precisely the 



conductance of a Hubbard chain of finite size Nc connected to noninteracting leads. |21l. |22 
the low-temperature Fermi-liquid regime of which had been studied before with a pertur- 
bation theory. It was confirmed that at half-filling the conductance through the four 
dots, Nq = 4, decreases exponentially with increasing U reflecting a Mott-Hubbard type 

even size, for odd Nc the Kondo 
resonance emerges at the Fermi level, and it contributes to the Unitary-limit conductance 
g = 2eyh. 

In quantum dots, the on-site potential is a tunable parameter, and it corresponds 
to the gate vol tag e. We have reported the dependence of the conductance of the triple 
quantum dots.[22| The purpose of the work is to study how the conductance through the 
four quantum-dot chain behaves as a function of the gate voltage. Owing to the Fermi-liquid 
properties at low temperatures, the conductance g and local charge A'ci in the quantum dots 
are determined by the two phase shifts Seven and ^odd, which are defined with respect to 
the even and odd [s and p) partial waves, .e.peet.ve,.Q We calculate tKe.e .J. pha.e 
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shifts from the fixed-point eigenvalues of NRG. Furthermore, from these phase shifts, we 
can also deduce a parallel conductance through the quantum dots connected transversely to 
four noninteracting leads. 

The results of the dc conductance show the typical Kondo plateaus of the Unitary limit 
g ~ 2e^/ h when the total number of electrons N^i in the four dots is odd. On the other hand, 
the conductance shows wide minima for even A'ci — 2, 4, 6. The local charge A'ci shows a 
quantized staircase behavior as a function of e^. Among the conductance minima, the one 
at half- filling A'ci ^ 4 is the widest, and it can be regarded as a kind of the Mott-Hubbard 
insulating behavior. The feature of the conductance near half-filling is quite different whether 
the number of quantum dots is even or odd, namely the Kondo plateau appears for odd Nc- 

This paper is organized as follows. In Sec. IHl we describe the outline of the formulation 
to deduce the conductance from the fixed-point energy levels of the NRG. In Sec. IIIIl we 
show the NRG results. In Sec. IIVI discussions are given. 



II. MODEL AND FORMULATION 



In this section, we describe briefly the relation between the phase shifts and fixed-point 



Hamiltonian of NRG.|22| As far as the zero-temperature value is concerned, the conductance 
can be calculated more accurately with this formulation than to calculate it directly from 
the current-current correlation function. 

We start with a Hubbard chain of a finite size as illustrated in Fig.[T] It is connected 
to two non-interacting leads on the left(L) and right (i?) via the tunneling matrix elements 

and f^, respectively. The Hamiltonian is given by 

H = He + Ti.Q + Timix + Hleaxi (l) 
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FIG. 1: Schematic picture of a series connection. 
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where di^ annihilates an electron with spin a at site i in the Hubbard chain, which is char- 
acterized by the nearest-neighbor hopping matrix element t, onsite energy e^, and Coulomb 
interaction U. In the lead on the left or right (z/ = L, R), the operator c]^^^ creates an 
conduction electron with energy eku, the wavefunction for which is denoted by (pkui'f')- The 
linear combinations of the conduction electrons t/;^^ and mixed with the electrons in 
the interacting sites at i = 1 and Nc, respectively, where ipua — Ylik^kua't'kuiju) and 
is the position at the interface in the lead v. We assume that the hybridization strength 
r,^ = T[v1 \(l>kv{rv)^ <^(^ — ^kv) is a constant independent of the frequency and take the 
origin of the energy to be /i = 0. 
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A. Ground-state properties and phase shifts 



In the following, we assume that the system has an inversion symmetry, vl = vp> (= v) 
and = (= r). Then, it is convenient to introduce the orbitals which have even and 
odd parities 

- ^2 ' ^ ' 

where j = 1,2,..., Nc/2 for even Nq- For odd Nc, there exists one extra unpaired orbital 
a(7Vc+i)/2,cr = c?(Arc+i)/2,(T in addition to the pairs labeled as j = 1, 2, . . . , {Nc — l)/2. Due to 
the inversion symmetry, at zero temperature T = and Fermi energy u = 0, the retarded 
Green's functions for ai^o- and 61,0- can be written in the forms 

«l,<xL=0 = f — = f / 2 11 ' 

J- l^even + « J- a/ + 1 

Namely, each of these two retarded Green's functions is determined by a single parameter, 
Keven or Kodd, which coutaius all effects of the scatterings and interactions. The phase shifts 
Seven and Sodd, Correspond to the angle of these two Green's functions in the complex plane. 



The two phase shifts determine the ground-state properties of the series connection illus- 
trated in Fig. ^[2^ Specifically, using the Kubo formalism and Friedel sum rule, the dc con- 
ductance (y'series and total number of electrons in the Hubbard chain = J2i^=i Ylia^Aadia) 



can be expressed in the forms, 



22| 

2e2 



5'scrics = ^ Sin^ (^5cvcn - ^odd j , (10) 
^el = - [Scvcn + ^odd) • (H) 

These formulas are justified if the imaginary of the Nq x Nc matrix version of the retarded 
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self-energy due to the interaction vanishes, Im S'^(O) = 0, at T = and = 0. 

Furthermore, from these two phase shifts defined with respect to the series connection, 
one can deduce the parallel conductance S'paraiici for the current flowing in the horizontal 
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FIG. 2: Schematic picture of a parallel connection. 



direction in the geometry illustrated in Fig. El 

5'parallcl — I Sm Op 



+ sin^ 5odd 



(12) 



This expression is equivalent to the multi-channel version of the Landauer formula. |2(t| and 



it can also be derived from the Kubo formula. [2^ |22| Since two conducting channels can 
contribute to the total current, the maximum possible value of (y'paraiici is 4e^//i. Note that 
the tunneling matrix elements between the leads and interacting part are taken to be w/ \f2 
as shown in Fig. |21 Because of the inversion symmetry between left lead 1 (2) and right 
lead 1 (2), an odd combination of the states from these two leads is separated, and the 
interacting site for i = \ [i = Nc) is coupled to an even combination by the matrix element 
V. Therefore, the ground state properties for the parallel geometry are determined by the 
phase shifts Seven and 6odd that are defined with respect to the series connection in eqs. (jHl) 
and (jni). Specifically, N^i for parallel connection is equal to that for the series connection, 
and is given by eq. (fTT|. 
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B. Fixed-point Hamiltonian and phase shifts 



The two phase shifts Seven and Sodd can be deduced from the fixed-point eigenvalues NRG. 
In the NRG approach, a sequence of the Hamiltonian if^r is introduced, by carrying out the 



logarithmic discretization for the conduction band with the half- width D, 
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171 . llSj l as 



Hn = A(^-^)/^ {K + n''c + + ) , (13) 

^^mix = ^ 1] ( flhJl^a + 4,<x/o,La) 

+ fU^Na,. + dkjo,R. ) ' (14) 

(7 

u=L,R (7 n=0 

^ ^ fn+l,ua fn,U(7 ~^ fn,va fn+l,ua ^ ) (l^) 

where v = ^J2DY Ap^jix , = ^ jz\jx log ^ 5 

l_l/A"+i 

= - , = . (16) 

v/1 - 1/A2"+Vl - 1^2^^+^ 
As in the case of the single Anderson (or Kondo) impurity, Q, 0, the low-lying eigenval- 
ues of if AT for the finite Hubbard chain converge, for large A^, to the fixed-point values which 
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have one-to-one correspondence to the free quasi-particles of a local Fermi liquid 
The fixed-point Hamiltonian describing the free quasi-particles can be written in the form 

= A^^-^^/^ {h^' + ^mix + <]) , (17) 

where = + ^f^=iReE^(0) dl^dj^. The many-body corrections enter through the 
real part of self-energy at T = 0, = 0. With NRG, one can calculate accurately the quasi- 
particle energies e* with parity 7 (= "even" or "odd") from the fixed-point eigenvalues of 
if AT. Then, the parameters and Kodd defined in eqs. (jS)) and © can be deduced from 
the quasi-particle energy £* .|22i| 

K, = {v'/TD) lim i^A(^-i)/2g^(e;) . (18) 

Here, gj^^^uj) = ^^=q |v5m(0)P/(co' — e^) is the Green's function for one of the isolated leads 
described in i^i^f]. The eigenvalue and eigenfunction of for one isolated lead are denoted by 
em and ipm{n), respectively, for < n < A^. 



Nc = 4: U/{2TTt) = r/t = 0.12 




FIG. 3: The conductance fcries calculated using NRG method (circle dots) compared with the 
exact result(solid line). Here, U/{2-nt) = and T/t = 0.12. For NRG, we use A = 6.0 and 
t/D = 0.1. 

III. NRG RESULTS FOR THE GROUND-STATE PROPERTIES 

In this section, we apply the formulation described in the above to the four-site Hubbard 
chain, Nc = 4, attached to the reservoirs. In each step of the iteration, the Hilbert space for 
-ff/vT+i is constructed from that for H]^ and extra 16 states with respect to the two orbitals 
fn+i^Ra fN+i,La- liavc kept the lowest 1716 eigenstates for constructing the Hilbert 
space for the next step. With this procedure the discretized Hamiltonian Hn, which consists 
of 2(A^ + 1)+A'^c' sites, can be diagonalized exactly up to = for Nc = 4. We have checked 
that the conductance of the the noninteracting case U = are reproduced sufficiently well 
with this procedure also for a rather large value of the discretization parameter A = 6.0, 
as shown in Fig. IHl We have also confirmed that the numerical results for the fixed-point 
eigenvalues of can be mapped onto the energy spectrum of the free quasi-particles in all 
parameter sets we have examined. It justifies the assumption of the local Fermi liquid we 
have made in the previous section. 

A. Series connection 

We now show the dependence of the conductance of the quantum dots in the series 
connection. As already mentioned, the onsite potential is a tunable parameter that is 
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Nc = 4: u/{2Tit) = 1.0, r/t = 0.12 




FIG. 4: (a) The conductance ^series, (b) local charge A'ci, phase shifts 2(5cvcn/7r and 2(^odd/7r as 
functions of ed/U for U/{2TTt) = 1.0, T/t = 0.12, t/D = 0.1, and A = 6.0. The dashed vertical lines 
in (b) correspond to the values of e^, at which A'^ei jumps in the limit of F ^ 0. 



controlled by the gate voltage in real quantum dots systems. In Fig. HI (a) the conductance 
fi'series' 0^) ^otsl chaigc A^ci in the four-site Hubbard chain, the phase shifts Seven and 5odd 
are plotted as functions of e^/U, where the parameters are chosen to be U/{2Tit) = 1.0, 
r/t = 0.12, t/D = 0.1 and A = 6.0. The vertical dashed hues in (b) correspond to 
the values of at which N^i jumps discontinuously in the limit of F ^ 0. The Kondo 
plateaus of the Unitary hmit (7sci.ios — 2e^//i are seen in the regions corresponding to the odd 
occupations A'ci ^ 1, 3, 5, 7, where the value of the phase shifts is almost locked such that 
Seven — ^odd — 7r/2. On the other hand, for even occupancies A'ei — 2, 4, 6, the conductance 
is suppressed substantially, and the phase shifts are locked at 5cvcn — ^odd — or vr. The 
conductance valley at half-filling becomes wider and deeper than the other valleys. Although 
the size of chain is not so large, i.e., Nc = 4, the feature of the valley at half-filling can be 
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FIG. 5: (a) The conductance ^series, (b) local charge A'ci, phase shifts 2(5cvcn/7r and 2(^odd/7r as 
functions of ed/U. Here, U/{2TTt) = 0.2, T/t = 0.12, t/D = 0.1, and A = 6.0. The dashed vertical 
lines in (b) correspond to the values of at which N^i jumps in the limit of T ^ 0. 

understood as a kind of the Mott-Hubbard insulating behavior. We will discuss later the 
behavior of the conductance near half-filling more quantitatively. In Fig. EJ the hybridization 
energy scale F is chosen to be much smaller than the hopping matrix element t. Therefore, 
the local charge TVei shows a clear stair case step behavior. In the limit of F — * 0, the local 
charge jumps discontinuously at the values of corresponding to the dashed lines. Note that 
unless the Coulomb interaction, the phase shifts do not show the 7r/2 steps corresponding to 
the plateaus for odd A'ei = 1, 3, 5, and 7. When the onsite potential ed/U decreases from 0.4 
to 0, the phase shift for the even partial wave 5evon increases from to vr via two successive 
7r/2 steps, while that for the odd partial wave 5odd remains almost unchanged at 0. Then, 
in the region of —0.5 < ed/U < 0.0, 6odd shows one 7r/2 step, while the even part remains 
to be Seven — Similar features are seen also for ed/U < —0.5. In fact, the two regions 



10 



-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 



FIG. 6: The conductance near half-filling for several values of the Coulomb interaction U/{2TTt) = 
0.2, 0.5, 1.0 and 1.5. Here, V/t = 0.12, t/D = 0.1, and A = 6.0. 

ed/U < —0.5 and e^/f/ > —0.5 are relating each other via an electron-hole transformation. 
In the case of the linear chain of the quantum dots, the resonance states of the even parity 
for the even partial wave and that of the odd one cross the Fermi level alternatively. Thus, 
the 7r/2 steps took place for even and odd phase shifts alternatively. 

We next examine the ground-state properties for small U . The results of the conductance 
and local charge for f//(27rt) = 0.2 are shown in Fig. where T/t = 0.12 and t/D = 0.1. 
The peaks of the conductance becomes very narrow compared to those for larger U shown 
in Fig. m Nevertheless, the shape of the peaks seen deviates from a simple Lorentzian shape 
as a result of the many-body effects due to the small but finite U. The conductance valley 
at half-filling N^i ^ 4 is slightly wider than the other ones at N^^i ~ 2 and 6. However, the 
difference is not so large. In Fig. El (b), the steps for the local charge are seen clearly for 
even occupancies ^ 2,4, and 6. However, for odd occupancies, only a weak structure 
can be recognized just for N^i — 1 and 7. Therefore, when a resonance state passes through 
the Fermi level, two electrons occupy the state almost simultaneously. Correspondingly, the 
phase shifts 5cvcn and 5odd do not show clear 7r/2 steps, although the weak structure can be 
seen for N^i — 1 and 7. 

In order to see how the conductance valley at half-filling evolves with the repulsive in- 
teraction, the conductance is plotted for several values of U/(27rt) = 0.2,0.5, 1.0 and 1.5 in 
Fig. El where the other parameters are chosen to be T/t = 0.12, t/D = 0.1, and A = 6.0. 
As U increases, the valley becomes wider, and the edge at the both side becomes sharper. 
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Furthermore, the minimum value of the conductance at the electron-hole symmetric case 
ed = —U/2 decreases exponentially with increasing U, which has been confirmed precisely 
in a previous work. 2]J These results imply that for even Nq the conductance valley at half- 
filling approaches to the insulating gap, even though the size of the chain examined here is 
just Nc = 4. However, because Nc is finite, the low-lying energy states of the whole system 
including the noninteracting leads are described by the Fermi-liquid fixed point. Therefore, 
the ground state is not an exact insulating state, but a highly renormalized metallic state 
with heavy quasi-particles. For large even Nc, however, the conductance is expected to 
show an exponential dependence g oc e~^^^^, where ^ ~ hvp/Agi^p is a correlation length 
determined by the Hubbard gap Agap and Fermi velocity vp. 

So far, the hybridization energy scale F is taken to be much smaller than the hopping 
matrix element t between the dots, F = 0.12t. In order to study the ground-state properties 
for larger F, we have carried out the calculations taking the hybridization to be F/t = 0.5. 
The results are show in Fig. where f//(27rt) = 1.0, t/D = 0.02 and A = 6.0. The 
local charge A'ei, which shows a staircase behavior for small F, becomes a gentle slope in 
Fig. m (b). This is because the large hybridization makes the resonance peaks broad, and 
it reduces effectively the correlation effects. Therefore, the peaks of conductance become 
round. Furthermore, the conductance valleys especially the ones corresponding to A^^ei — 2 
and 6 become shallow. The valley at half-filling is still deeper than the other two valleys 
because of the correlations that lead the Mott-Hubbard behavior in the thermodynamic 
limit. 



B. Parallel connection 

We have also calculated the parallel conductance Sfparaiici using eq. (fT^ with the phase 
shifts obtained for the series connection. The results are shown in Fig.|Hlas functions of e^/t/ 
for U/{27it) = 1.0. In the figures, both the series and parallel conductances are plotted for 
comparison, and for the hybridization two different values are examined, (a) F/t = 0.12 and 
(b) r/t = 0.5. For small F, the parallel and series conductances are almost coincide with 
each other as seen in Fig. |Hl(a). The difference becomes larger, however, for large F, as seen 
in the lower panel (b). The valleys of the conductance are deeper for the parallel connection 
than that of the series connection. This seems to be caused by the fact the even and odd 
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FIG. 7: (a) The conductance fericsi (b)local charge A'ci, phase shifts 25cvcn/T^ and 2(5odd/7r as 
functions of e^/f/. Here, U/{2TTt) = 1.0, T/t = 0.5, t/D = 0.02, and A = 6.0. The dashed vertical 
hues in (b) correspond to the values of at which N^i jumps in the hmit of T — > 0. 

parts of the partial waves contribute to the parallel conduction separately through eq. (|T2|) . 
The parallel conductance for large F, shown in Fig. |Hl(b), is slightly larger than 2e^/h at the 
first and fourth peaks. As mentioned in the previous section, the upper bound of Qparaiici is 
2e^/h X 2, because two conducting channels can contribute to the current in the geometry 
shown in Fig. |21 However, in the present case, the parallel conductance does not reach the 
upper bound. Similar feature has also been confirmed to be seen for the triple dots in a 
previous work. As mentioned, in the case of the linear chain of the quantum dots, the 
resonance state of the even parity and that of the odd parity passes through the Fermi level 
alternatively when increases (or decreases). Therefore, only a single resonance state can 
contribute to the tunneling current for given e^, so that the peak value of the conductance for 
the linear chain is bounded at (y'paraiici — 2e^//i despite having the two conducting channels. 
To reach the Unitary-limit value 4e^//i for two conducting channels, the quantum dots are 
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iVc = 4 [//(27ri) = 1.0 




FIG. 8: The series (solid line) and parallel (dotted line) conductances for U/{2Trt) = 1.0 as 
functions of Cd/U. Here, (a) T/t = 0.12, t/D = 0.1, and (b) T/t = 0.5, t/D = 0.02. 

required to have a degeneracy in the discrete energy levels. It is not the case of the linear 
chain considered in the present work. 

IV. DISCUSSIONS 



As discussed in Sec. Illlj a kind of the Mott-Hubbard insulating behavior is seen at 
half-filling in the conductance through the one- dimensional array consisting of the even 
number of quantum dots Nq- The situation is quite different for the array with odd A^^;, 
because the Kondo resonance contributes to the current as shown in Fig. El where the 
series conductance and total charge Nd in triple dots {Nc = 3) are plotted as functions of 



ed/U. 



22l | Near half-filling, the conductance shows a plateau of the height 



2e^/h in 



stead of a wide gap Agap as seen in Fig. IHl for even Nc- The plateau becomes wider with 
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Nc = S T/t 
1 




FIG. 9: NRG results for the conductance feries and local charge A'ci for triple dots Nc = 3 as 
functions of ea/U. Here, U/{2Tit) = 1.0, T/t = 0.12, t/D = 0.1, and A = 6.0. The dashed vertical 
lines correspond to the values of at which A'ci jumps in the limit of T — > 0. 



increasing [/.j22] However, the Kondo temperature Tk decreases with increasing ?7,|2]J and 
thus the conductance plateau can be observed only at low temperatures T < T^- If the 
single-particle spectral function is calculated at half-filling as a function of the frequency 
cj, the narrow Kondo peak will be seen at < |ci;| < Tk- Then, there will be a pseudo- 
gap region at Tk ^ \^\ < Agap, where the spectral weight almost vanishes. The width 
of the Kondo plateau at half-filling will be almost the same as the excitation gap 2Agap 
seen in the spectral function. The Kondo temperature Tk for half-filling decreases with 
increasing Nc, and finally the resonant state vanishes as Tk — in the limit of Nc oo- 
Therefore, for large Nc, the conductance plateau emerges only at very low temperatures. 
If the thermodynamic limit A'^c' — > oo is taken first at finite temperature T ^ and then 
the limit T — > is taken next, the conductance will vanish because of the Mott-Hubbard 
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EOM results 
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U/(27it)=1.0, r/t=o, 
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FIG. 10: (Color online) EOM (equation of motion) results for the series conductance for the 
quantum-dot array of the size (a) Nc = 3 and (b) Nc = 4 obtained from eq. H20|) . where U / (27rt) = 
1.0 and T /t = 0.12. The results describe a typical feature of a Coulomb oscillation at high 
temperatures Tk <^ T <^ U . The dashed lines show the local charge A'ci in the isolated molecule 
for r = 0. 



behavior. 

So far, we have concentrated on the transport properties at zero temperature. For com- 
parison, we now consider the conductance at high temperatures. For quahtatively under- 
standing of the physics at high-energy energy scale, the equation motion (EOM) method 
can be used.^^ It is equivalent basically to the Hubbard I approximation, Q] and in the 
present case we start with a molecule limit, where F = and the quantum-dot array can be 
regarded like an artificial molecule described by the Hamiltonian 7i^ + 7i^ defined in eqs. (0) 
and Q. With this method, the full Green's function G^^'~'^\uj) is obtained by substituting 
the self-energy defined with respect to the molecule S^™°'''(u;), which is calculated exactly, 
into the Nc x Nq matrix Dyson equation, 

|^(EOM)^^^|-i = |g(°)m}~' - . (19) 



Here, G^'^\uj) = {G^^^\uj)} is the noninteracting Green's function defined with respect to the 
whole system including the two leads, and is determined by the Hamiltonian + T^mix + 
TYicad- Then, a typical value of the conductance at high-temperatures Tk -C T -C [/ is 
estimated by using an approximate formula following Kawabata,!^] 



e 



2 



— X 4F2 
h 



<r(0) , (20) 



2 



where a factor 1/2 has been introduced phenomenologically taking into account the fact that 
at high-energies the resonance tunneling using the Hubbard band occurs not simultaneously 
for the up and down spin components. In Fig. the conductance obtained from eq. ((201) 
plotted for (a) Nc = 3 and (b) Nc = 4, where the same parameter value is chosen for U and 
F as those in Figs. El and lU At high temperatures, the conductance almost vanishes both for 
even and odd A'ci. The peaks of the Coulomb oscillation appear at the values of e^, where 
iVei jumps discontinuously by an addition of one electron. If the temperature decreases, the 
bottom of the valleys at odd Nc,\ will rise, and it develops at T < Tx to the Kondo plateaus 
of the Unitary limit which we have described in the previous section. Particularly, the wide 
gap for the triple dots near half-filling seen in Fig. ^1 (a) evolves to the broad plateau in 
Fig.ini 

In summary, we have studied the ground-state properties of an array consisting of four 
quantum dots based on a Hubbard chain attached to two non-interacting leads. Using NRG 
approach, we have deduced the phase shifts, by which the conductance and local charge 

17 



away from half-filling can be determined, from the low-energy eigenvalues near the Fermi- 
liquid fixed-point. We have also discussed the parallel conductance of the quantum-dot array 
connected transversely to four leads. Our formulation to calculate the two phase shifts for 
even and odd partial waves is quite general, and will be applied to various quantum-dot 
systems in the Kondo regime. 
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